Gravitational waves and dragging effects 

Jifi Bicak^'^* Joseph Katz^'^'l" Donald Lynden-BelF* 

^Institute of Theoretical Physics, Charles University, 180 00 Prague 8, Czech Republic 

"^The Racah Institute of Physics, Civat Ram, 91904 Jerusalem, Israel 

OO ^Institute of Astronomy, Madingley Road, Cambridge CBS OH A, United Kingdom 

^ July 19, 2008 



I Abstract 

Linear and rotational dragging effects of gravitational waves on local inertial frames are 
studied in purely vacuum spacetimes. First the linear dragging caused by a simple cylin- 
drical pulse is investigated. Surprisingly strong transversal effects of the pulse are exhib- 

I ited. The angular momentum in cylindrically symmetric spacetimes is then defined and 

> 

^vj confronted with some results in literature. In the main part, the general procedure is 

^^ developed for studying weak gravitational waves with translational but not axial symme- 

try which can carry angular momentum. After a suitable averaging the rotation of local 
l> 

^^ inertial frames due to such rotating waves can be calculated explicitly and illustrated 

00 
^^ graphically. This is done in detail in the accompanying paper. Finally, the rotational 

. 1^ dragging is given for strong cylindrical waves interacting with a rotating cosmic string 

5_j with a small angular momentum. 
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1 Introduction 

Bets have been won by those who prognosticated that gravitational waves will not be pos- 
itively detected in an earth-bound laboratory until 1999, 2001,..., 2008. But it would be 
hard to find a reader of Classical and Quantum Gravity who does not believe that eventually 
gravitational waves will be detected through their effects on test laboratory masses. From 
purely theoretical perspectives, from the year 2000 at least six papers appeared in this jour- 
nal in which dragging effects even on gyroscopes were studied: for example, the precession 
of a small test gyroscope induced by a weak plane gravitational wave [1] , [2] , the behaviour 
of a gyroscope located at smooth and polyhomogeneous null infinity [3], [4], the motion of 
spinning test particles in plane gravitational waves [5], [6]. 

However, a truly "Machian" problem associated with gravitational waves is not concerned 
with the behaviour of a gyroscope immersed directly in a gravitational wave. A collapsing, 
slowly rotating spherical massive shell drags into rotation inertial frames inside its flat inte- 
rior [7], [8]. In linearly perturbed Friedmann-Robertson- Walker universes, the local inertial 
frames can be seen to be determined instantaneously via the perturbed Einstein field equa- 
tions even from the distributions of energy, momentum and angular momentum beyond a 
cosmological horizon [9], [10]. Correspondingly, it is important to demonstrate whether en- 
ergy and angular momentum of gravitational waves in purely vacuum spacetimes can cause 
the local inertial frames to rotate. If, indeed, this is so, it is of interest to see what such 
effects look like explicitly. 

Thanks to progress in "analytic gluing techniques" in geometric analysis and differential 
geometry over the past 25 years, and their recent applications to General Relativity, the 
question of the existence of the effect of the angular momentum of gravitational waves on 
the rotation of inertial frames can now be answered with full rigour. In our context the 
important physical statement following from these developments, in particular from the work 
of Corvino [11] and Corvino and Schoen [12], is that asymptotically flat initial vacuum data, 
i.e., data which represent pure gravitational waves can be deformed so that in the far region 
they are exactly Kerrian, i.e., stationary. Since the Kerr metric is one of the best prototypes 
to demonstrate the dragging of local inertial frames , there is now clear theoretical evidence 
available that angular momentum due to gravitational waves causes local inertial frames to 
rotate with respect to the local inertial frames at infinity. Moreover, as in the cosmological 



perturbation theory, or in the examples with rotating shells, this effect is global and instan- 
taneous (all the construction occurs on a given spatial hypersurface) . From the perspective 
of Mach's principle, one would prefer to demonstrate the rotation of local inertial frames in 
a dosed universe filled only with gravitational waves since then boundary conditions (like 
asymptotic flateness) play no role. Such a study lies in the future. However one would also 
like to see the dragging effects caused by the "rotating waves" explicitly. To construct explicit 
examples, we turn to the waves with a symmetry along Oz in cylindrical polar coordinates. 

In the next section the metric with a translational symmetry is introduced in suitable 
coordinates and the vacuum field equations in four dimensions are rewritten as Einstein's 
equations in three dimensions with a scalar field as a source, following [13]. A rotational 
symmetry is then assumed in addition, and the equations for standard Einstein-Rosen waves 
are written down, including one of their simplest explicit solution, the Bonnor/Weber- Wheeler 
time symmetric pulse. Section 3 is devoted to the effects this pulse exerts on test particles. 
Specifically, we calculate the force acting on a particle which is at rest on a cylinder with 
either fixed proper circumference or a fixed proper radial distance from the axis. Interestingly, 
it can happen that no physical force can keep particles at fixed circumferences in the "strong 
region" of the pulse. However, the particle can be forced to stay fixed at a given proper 
distance from the axis or at a given radial coordinate (defined geometrically in Einstein- 
Rosen spacetimes). The "linear stretching" due to the cylindrical wave in the directions 
transverse to its propagation thus appears to be stronger than the "linear dragging" in the 
radial direction. 

In the following sections we turn to the "rotational dragging". First, however, we ana- 
lyze the angular momentum in cylindrically symmetric spacetimes (Section 4). Those space- 
times posses the rotational Killing vector which can be normalized at the axis of symmetry and 
used to define the total angular momentum density, contained in a 2-dimensional cylinder of 
a unit length, by adopting the Komar expression. The Komar-type expression yields mean- 
ingful results even if a rotating cosmic string is present along the axis. We show that among 
cylindrical metrics which appeared recently in the literature the angular momentum is non- 
vanishing only in those cases which represent cylindrical gravitational waves interacting with 
a rotating cosmic string. 

An interplay between a strong gravitational wave and a rotating string is of interest but a 
physically "clean" explicit demonstration of dragging effects due to gravitational waves should 



not involve any matter like a string. A simple procedure for studying gravitational waves with 
a net angular momentum is presented in Section 5. We assume the waves to have the trans- 
lational symmetry but not axial symmetry. Then it is convenient to make dimensional re- 
duction after which the vacuum problem in four dimensions becomes that of a scalar field 
in 3-dimensional spacetime [13]. The scalar field satisfies the wave equation in flat space- 
time provided that it is weak. We can construct "(/^-dependent" , rotating solutions of this 
flat-space wave equation and obtain the effective energy-momentum tensor - which is of the 
second order - on the right-hand side of the Einstein field equations in three dimensions. The 
resulting metric deviates from 3-dimensional Minkowski metric by terms of second order. 
The rotation of inertial frames on the axis cannot depend on Lp. We are therefore primarily 
interested in the behaviour with radius of the (^-averaged dragging rate. We discuss the struc- 
ture of all the field equations, their consistency and suitable boundary conditions within this 
approach. The same function in the metric which approaches a constant at spatial infinity 
and represents the total energy in the case of the Einstein-Rosen waves, also yields the energy 
of rotating waves. The explicit solution of these equations based on a non-axially symmetric 
generalization of a weak Bonnor/Weber- Wheeler pulse, are constructed in the accompanying 
paper [14] . Therein the dragging effects of the resulting gravitational waves on local inertial 
frames are analyzed in detail and graphically illustrated. In Appendix A, the Ricci tensor 
for the averaged rotating gravitational waves in 2-|-l dimensions is given. 

The rotational dragging in the spacetime of the exact (in general strong) Bonnor/Weber- 
Wheeler cylindrical pulse perturbed by a rotating cosmic string with a small angular momen- 
tum is considered in the last Section 6. 

2 The metric and field equations for axially rotating gravita- 
tional waves 

At the core of this work and the following paper [14] are source-free metrics with at least one 

hypersurface orthogonal Killing vector of spacelike translations C,. In this case coordinates 

exist ^ 

{x''} = {x",a;3 = z} , where C = {0,0,0,1}, (2.1) 

^Greek indices A, ^.,v,p,- ■ ■ run over the four spacetime coordinate labels 0, 1, 2, 3; Latin indices a,b,c,--- 

run over the time and two spatial coordinate labels 0, 1,2. The metric g^,/ has signature H and g is its 

determinant. Covariant derivatives are indicated by a D, partial derivatives by a 9 and covariant derivatives in 
a three subspace by Va- The permutation symbol in 4 dimensions is e/j^^pa- with eoi23 ~ 1, Vf^vpcr = x/^Em^po-- 



and the metric can be written in the form introduced by Ashtekar, Bicak and Schmidt [15]: 

ds2 = e-'^'^gabdx'^dx^ - e^'^dz^, (2.2) 

ip and Qab are functions of x'^ only. The source-free Einstein's equations R^y = take the 
following interesting form in terms of the Ricci tensor IZah of the 3-space da"^ = gabdx'^dx^: 



= ^ ^Zab = 2^aiJ^b^l^, (2.3) 

R33 = => 5"'VaVfeV = 0. (2.4) 

These equations can be interpreted as Einstein's equations in three dimensions with a scalar 
field J- = ip/^/AirG as a source. 

The Einstein-Rosen metrics [16] represent a particular class of metrics of the form (2.2) 
which admits an additional hypersurface orthogonal Killing field of rotations t}. In this case, 
coordinates are generally chosen such that 17= {0,0,1,0} and are denoted by x^ = t,x^ = 
p, x^ = ip. The metric then reads 

ds" = e-''^ [e^^ [dt^ - dp") - p^d^^] - e^'^dz^ (2.5) 

ij) and 7 are functions of t and p. Notice that the radial variable p is a geometric object 
[17]. In this case the non-trivial equations (2.3) and (2.4) can be written as follows (in the 
formulas, d^X = dtX = X and diX = dpX = X' for any X{t, p)): 

7^oo = 7" + -7-7 = 2V'^ (2.6) 

P 

7^ll = -7" + -7' + 7 = 2V'', (2.7) 

P 

7^ol = -7 = 2^^', (2.8) 

P 

and 

< + -^' -^ = {). (2.9) 

P 

The last equation is that of cylindrical waves in cylindrical coordinates in flat space. Given 
a solution ^ satisfying appropriate boundary conditions, the function 7 is entirely defined by 
the following combinations of equations (2.6) - (2.8): 

7' = p J'^' 2 + ^A and 7 = 2pV'V''- (2-10) 

The wave equation (2.9) guarantees that the integrability conditions of the last two equa- 
tions are satisfied. 



Of particular interest are the Bonnor [18] and Weber- Wheeler [19] time symmetric in- 
coming and outgoing waves which are smooth and finite everywhere at all time. They have 
been also discussed in some detail in Weber's book [20]. That metric is not flat at spatial 
infinity. Its asymptotic behavior received special attention in [15]. The explicit forms of ip 
and 7 given below are taken from Ashtekar, Bicak and Schmidt [13]; we shall write them 
with non-dimensional variables: 

, V2C 
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(2.11) 
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C and a are two constants of integration: a is a measure of the width and C a measure of 
the amplitude of the waves^ [20] , and [19] : 
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1 illustrates ^(t, p) and 7(t, p) as functions of p for various values of i. Notice that 
for p ^ CO. 




Figure 1: The plots of 'i/'(t, p) for |t| — 0, 1, 2, 5, 9 from the left to the right, and of 7(t, p) as functions 
of p for |i| = 0, 2, 5, 9 from the left to the right. 



^We shall refer to them simply as "the width" and "the amplitude" 



3 The linear dragging in a Bonnor and Weber- Wheeler wave 

One way to calculate the acceleration due to gravity on earth is to calculate the force per 
unit mass necessary to oppose it and prevent the body from falling. To a relativist that force 
has the spatial components of the 4-acceleration necessary to keep the body at rest. Here 
we evaluate the force necessary to keep a particle in various geometrically defined positions 
subject to a Bonnor/Weber- Wheeler pulse. Surprisingly we find it sometimes impossible to 
keep the particle on a cylinder of a fixed proper circumference length. However, it is possible 
to keep the particle at a fixed proper distance from the axis or at fixed p. These worldlines 
are always timelike. 

Let v^^ = dx^/dr, QfivV^v'^ = 1 represent the 4-velocity of a particle with proper time r 
which is not falling freely. Notice that the momenta of a freely falling particle in the 99 or z 
directions are constants of motion because of the global symmetries of spacetime. Thus any 
acceleration in the ip 01 z directions can only be due to non-gravitational forces. Although the 
wave can stretch significantly distances between points in transverse directions, gravitational 
forces act only in the p direction. The interesting cases here are thus those of particles forced 
to stay or move in the p direction. Their 4-velocity has two non-zero components v^ and v^ . 
Following (2.5), 

g^^v^v'' = e2(^"'^) [{vy - [v^] = 1, (3.1) 

while fo = e^^'^^'^'v^ and vi = —e'^^'^~''^'v^. So equation (3.1) can also be written as 

(^o)2-(t;i)2 = e2(''-^). (3.2) 

The momenta per unit mass V2 = v-^ = 0. The 4-acceleration is a^ = Dv^/Dt = dv^/dr + 
2 VpVo-d^g^'^. It has two components, 

ao = -} ii - Ip) and ai = — (7' -7/;'). (3.3) 

ar dT 

However, if the particle is forced to move along the worldline [t(r),/9(r)], so that v^ = 
w^[t(r),/9(r)], we have successively 

—— = vqv + vq V and — — = viv + vi v . (3.4) 

dr dr 

On the other hand the derivatives of (3.2) give respectively 

vqVq — vivi = {^ — 'il>)e'^~^ and vqVq — vivi = {^' — ijj')e^~'^ . (3-5) 



We may eliminate for instance vq and vi from the pair of equations (3.4) and (3.5). The new 
expressions for dv^/dr and dvi/dr may be substituted in (3.3) which becomes 

flQ = f {vq — vi) and ai = —v {vq — vi). (3.6) 

This makes it obvious that a^v^ = 0. 

The frame component of the force / per unit mass or the acceleration exerted on the 
particle is 

/ = ^/=5ira' = -V^g^ai = -e-^-'-^^ai. (3.7) 

The magnitude of the acceleration is 

yZ^^ = e-2(7-V')|^;g'_^^|. (3.8) 

Suppose now a particle sits at (p = (po and z = zq and its worldline is given by 

F{t, p) =C , C is constant. (3.9) 

If F' 7^ 0, we may use C instead of p as a spatial coordinate. Differentiating (3.9), we obtain 
dC = Fdt + F'dp, from which follows that dp = {l/F'){dC — Fdt). Substituting this dp into 
(2.5) we obtain 



ds 



2 ^ 



pi2 



F'^ -F^)dt^ + 2FdCdt-dC^ - e~'^'^ p^ dp^ - e^'^ dz"^ . (3.10) 



The particle in question is at rest in {t, C, p, z} coordinates. The proper time elapsing after 
an interval dt is thus dr = (e'^~'^ /\F'\) \/F''^ — F'^dt, and the particle can only be at rest if 

F''^-F^>0. (3.11) 

Otherwise C = F{t, p) represents a spacelike worldline. 

Assuming the inequality (3.11) is satisfied, we may calculate the force per unit mass that 
must act on a particle to be kept at fixed C. Since 

dF 

^ = Ft,o + F'v^ = e-^^^-'^^Fvo - F'vi) = 0, (3.12) 

dr 

equation (3.12) with (3.2) define vq and vi. One readily finds that \i vq > {) which is always 

allowed, 

^^-■^ 
Vq = epiF'H , vi = epiFH with epi = sign{F') , H = — =, (3.13) 

VF'^ - F2 



and the force per unit mass is given by (3.7) or 



/ = e^-V(V-i;i). 



(3.14) 



Now the particle will stay on a circumference of length I ii F = pe~^ = 1/2tt, or in 
dimensonless form, L = pe~^ = l/2iTa. Then, dpL = (1 — ptlj')e^^ and d^L = —pipe'"^, and 
condition (3.11) amounts to having (dpL)'^ — {d^L)^ > 0, or 



u{t,p) = {i-pi^r-{pi^y>o- 



(3.15) 



We are inter U(i,p) as a function of L{t,p). We consider first an amplitude 6 = 1. We shall 
see later what happens when b is not equal to 1. 

L{t,p = 0) = and L{t,p = cxd) = cx). L increases monotonically for any |t| less than 
about 6. For |t| > 6, L decreases in some interval say {pi,p2} which depends on \i\ and is 
always around p = \t\; see Figure 2. 




Figure 2: From left to right L as a function of p for \i\ — 1,6, and 15. 

Now consider U. U{t,p = 0) = 1 and U(t,p = oo) = 1. Thus U is positive for small 
and big L's. What happens for intermediate L's? Consider U as a function of L for various 
values of the dimensonless time t. For \t\ < 2, [/ > 0. U varies monotonically when p 
nears and oo. A typical example is given in Figure 3a where \t\ = 1. Notice that in 
the graphs of Figure 3 the U axis cuts the L axis where \t\ = p. This is done for reasons 
appearing later. Figure 3b shows how U{L) varies when |t| ~ 1.95; this is a limiting case. 
It becomes impossible to keep a particle on some of the cylinders at \t\ > 1.95. For instance 
for |t| = 6, when L has a horizontal inflexion, U < for 4.22 ^ L < 4.56 corresponding 
to an interval 4.83 ^ p ^ 6.32. For |t| > 6, U{L) has a loop because L does no more 
increase monotonically; see for instance Figure 3c in which |t| = 20. At that time it becomes 
impossible to keep a particle at 17.05 ^ L < 17 A corresponding to 18.1 ^ p < 20. The 
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Figure 3: C/ as a function of L. A particle can stay at a given L only if [/ > 0. See text for details. 

loop persists at any later (t > 6) or earlier (t < —6) times and the knot is always where 
[/ > while U is always negative for |t| = p. See for instance Figure 3d for |t| = 1000. 

What if the amplitude b is different from 1? For b ^ 1 the phenomena is essentially the 
same. Loops appear at bigger L's for b < 1 and at small L's for 6 > 1. It may be noted that 
the U axis, which cuts the L axis where |t| = p also cuts the loop in the region where U < 0. 
It is therefore interesting to evaluate the upper bound on b resulting from the condition (3.15) 
for \t\ = p. In this way we find that U is always negative when \t\ = p \i 



b> bM 



(l+V)^l + (l + 4p2 



2p2 



(3.16) 



Thus, the phenomena described for 6=1 exists at any amplitude however small since for 
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p ^ oo, 6m ^ 25/Vp^/2 ^ 0. 

If a particle is kept at a fixed proper distance A =const from the axis, then condition 
(3.11) becomes 

-e^-"^ <dt(r e^-^dp^ < e^--^. (3.17) 

The situation is even clearer, and much easier to manage numerically, if the particle is kept 
at fixed^ p. Since p is a spatial coordinate there is manifestly no condition like (3.11). In this 
case, v^ = 0, wq = e'^~^, and the dimensionless force is given by / = a/ = —dpe~^"'~'^>. In 
real value, / = — /. The particle can always be kept at fixed radial coordinate p. 

As an example, let / ~ 0.3, a characteristic value for \t\ > 20 and 6=1. Let / be 
barely bearable, say f = bg "^ b ■ 10'^ cm sec^^. Then, a = ^ ~ 3 • 10^ AU. The amplitude 
\/2C = ab is of the same order as a. The characteristic time tc during which the acceleration 
must be withstood is tc" ~ - ~ 2 • lO^sec ~ 3weeks. 

4 Angular momentuin of cylindrically symmetric spacetimes 

Here we consider spacetimes with two spacelike Killing fields C, and r]. ^ is in general hyper- 
surface orthogonal but rotational rj is in general not, so that the two-dimensional isometry 
group is not orthogonally transitive [21]. The spacetime may admit an additional Killing 
vector ^ which is timelike and may thus represent a stationary field of a rotating massive 
cylinder or an infinite straight rotating cosmic string. 

Although in the main part. Section 5, we require the axis of symmetry to be regular, 
in Section 6 we also discuss some aspects of cylindrical spacetimes containing an infinitely 
thin string in the "wire approximation". It is associated with specific conical singularities 
along the axis which, though being "quasi- regular" , do not strictly speaking, belong to the 
spacetime. An extensive literature on cylindrical waves interacting nonlinearly with cosmic 
strings is available^. The wave solutions in vacuo are considered by some authors as too 
restrictive. The reason is a generalization of Papapetrou's theorem from stationary axisym- 
metric vacuum spacetimes to the cylindrical case: if the spacetime contains at least part of 
a regular axis of cylindrical symmetry, the isometry group must be orthogonally transitive, 
so Killing orbits admit orthogonal surfaces [23], [24], [25]. The orthogonal transitivity thus 

^See [17] where it is shown that p is a geometrically well defined quantity. 
''See the review [22] and references below. 
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excludes the possibility of a global rotation. In other words in vacuo there can be no 'rotat- 
ing cylindrical waves'. With a material source present as in the case of the rigidly rotating 
dust cylinder [26], for example, the spacetime can, of course, be regular everywhere with 
a non-vanishing angular momentum per unit length. Bondi [27] studied general changes in 
time of such systems which can lead to radiation. As he noticed, the conservation of an- 
gular momentum occurs even if gravitational waves are emitted by the cylinder since the 
cylindrical symmetry of the waves precludes their carrying angular momentum. 

In all the cases mentioned, the spacetimes posses the rotational Killing vector rj with 
closed orbits. This can be normalized at the axis of symmetry, with corresponding modifica- 
tion in case of a string present along the axis. When the rotational Killing vector is available 
we can define a total angular momentum density contained inside a 2-dimensional cylinder C 
of 'unit coordinate height', generated by orbits of the isometry group, by adapting the Komar 
[28] expression^ to cylindrical symmetry: 

JiC) = ^ f Dl^r^-^ldS^^, where dS^^ = \ e^^p^dxP A dx'' , t]^ = y=^r?^ , k=^ 
^1^ Jc 

(4.1) 

Since 

r = Y^d, (Z)[^r}'^]) (4.2) 

is conserved, d^j^ = 0, in vacuo the integral (4.1) is independent of the choice of C; the integral 

over 'bottom' and 'top' of the solid cylinder formed by two cylinders C and C vanishes due 

to cylindrical symmetry. This corresponds to the well known fact that the Komar energy, 

defined as (4.1) with -q being replaced by the timelike Killing field ^ in stationary spacetimes, 

is independent of the choice of the closed 2-surface as it is moved continuously through a 

vacuum region [30] . Hence, we can define the total angular momentum per unit length in the 

C, direction by 

J = ^ I D^'^r^d^^,, (4.3) 

2k Jc 

where Coo denotes the cylinder of unit length at infinity (there the proper and coordinate 
lengths along the C, direction coincide). 

Let us express the angular momentum (4.3) ("density" will often be omitted) explicitly 
for some specific metrics. In all cases the coordinates adapted to the symmetry are those 
defined in Section 2 below equation (2.4). The 2-surface element 

dSoi = — (iSio = d(pdz. (4-4) 



■^For using Komar's expression in case of the angular momentum, see [29] and [30] p. 296. 
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Hence the angular momentum (4.3) reads 



J =- I D^^fi^^difdz. (4.5) 

K r 



The integrand does not depend on ^p and z. Integrating we thus get ® 

27r 



J 



i^[o^i] 



(4.6) 



Mashhoon et al [31] considered a class of "rotating gravitational waves" which should 
represent "radiation that propagates outward or inward and at the same time has non-trivial 
azimuthal motion" . Their metrics satisfying vacuum field equations have the form 

ds" = e^^-^^{dt^ - dp") - W'e-^^{d^ + udtf - e'^dz^; (4.7) 

the functions 7, ip, W and uj depend on t and p only (in [31] p = R and W = p). If u; = 0, 
the field equations imply that W can be chosen equal to p, and the metrics (4.7) become 
the well known metrics of Einstein-Rosen waves (e.g. [21], [15]). Calculation of the angular 
momentum (4.6) yields 

J = -W^e-'^^uj' . (4.8) 

In Section 6 we shall especially need the angular momentum for small lo. One of the 
field equations (cf. [31] or the Appendix where the Ricci tensor in the reduced 3-dimensional 
spacetime is given) then implies W = p (like in the case of Einstein-Rosen waves), and 
Lu = zuoo/p'^ + 0{p~^). The angular momentum then becomes 



K 



27r K 

—e-^^^woo or vjoo = -^Je^^-'. (4.9) 

K zvr 



Next, consider the case of infinitely thin cosmic strings. When they posses an angular 
momentum J, the corresponding metric can be written as (e.g. [32], [33]) 

ds'^ = {dt + ^Jdiff - {p)-^f'/^{dp'^ + p^dy^^) - dz\ (4.10) 

ZTT 

p being the mass per unit length. After putting {j))^'^^d~p = dp, the metric can be rewritten 
as 

^^2 ^ (^^ ^ ^Jd^f - dp" - (1 - ^) W - dz\ (4.11) 

27r 2iT 



®Let us remark that with z replaced by 61 G [0, tt] and {t, p, 9, (p} denoting the Boyer-Lindquist coordinates 
in the Kerr metric (e.g. [30]) with mass M and angular momentum per unit mass a, the expression (4.6) 
yields indeed J — Ma. 
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In both metrics ip G< 0,27r)^. Comparing (4.10) and (4.11) with the metric (4.7) and 
the angular momentum (4.9) for small cj, resp. J, it is seen immediately that the Komar 
expression gives indeed J appearing in the metric. 

Various authors [34], [35], [36], [37] constructed new families of exact solutions of the 
vacuum Einstein (and Einstein-Maxwell) equations for cylindrically symmetric non-stationary 
spacetimes which they interpreted as gravitational (and electromagnetic) cylindrical waves 
interacting with a rotating cosmic string. In their solutions, however, the cross-term dtdip is 
missing. Applying our expression (4.6) for the angular momentum to their metrics, we get 
J = 0. More recently, Manojlovic and Marugan [38] considered thoroughly cylindrical waves 
interacting with spinning cosmic strings within the Hamiltonian formalism. In their work the 
non-diagonal terms proportional to dtdip are admitted ~ their value of the angular momentum 
coincides with that following from (4.9). Let us also state that our expression (4.6), taken 
at fixed p, when applied to Bondi's cylindrical metric (equation (1) in [27]) given in the 
Weyl-Papapetrou form, yields precisely the Bondi expression for the angular momentum. 

To the end of this section let us mention one plausible feature of the angular momen- 
tum defined above. Imagine we transform the metric (4.7) into rotating axes by introducing 
Tp = ip + 0,t, in which 0, may depend on time. Then dp = dp + 0,(t)'dt, but uj' = uj' so 
that the total angular momentum (4.9) remains unchanged. It thus characterizes an intrinsic 
property of the field. 

5 Rotating waves in the symmetry reduced general relativity 

To construct rotating gravitational waves with a non-vanishing angular momentum we can 
keep the translational symmetry but we have to drop the assumption of axial symmetry. In 
general we would have to solve the coupled system of equations (2.3) and (2.4) with ^ and gah 
depending on all three coordinates x" = t,x^ = p, x^ = p. In 3-dimensional spacetimes the 
Riemann tensor 

T^ahcd = 2 \{Tl-a{c - z9a[cT^) 9d]b " (^6[c " 4^9b[c'^) 9d]a\ ■ (5-1) 

The Ricci tensor can be determined directly from the field iJj by equations (2.3). 



^The transformation '(p — (1 — ff )v accompanied with the change of time t = t + ^J'^ leads to the 
flat metric with the range ^ £ [0, (f — |^)27r), corresponding to conical geometry and the time coordinate 
'jumping' by kJ when the string is circumrotated. 
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Clearly to tackle such a problem is a formidable task. However, the above form of Ein- 
stein's equations suggests naturally the following approximation procedure. Assume the field 
tp and its derivatives to be small, V' = €^{x'^), e is a small dimensionless parameter. Then 
the field equations (2.3) show that TZab ~ C(e^) and, according to (5.1), the Riemann tensor 
can be calculated in terms of datp up to ©(e^). We may thus write the 3-dimensional metric 
in the form 

gab = Val> + ehal>{x') or g-^ = r^-^ - e^ ^^x') . (5.2) 

Therefore we can construct a rotating ("(/^-dependent") solution of the wave equation in flat 
space in cylindrical coordinates {t,p,ip} and still satisfy the field equations (2.4) in terms 
~0(e2). 

Next, we can write down the Ricci tensor on the left-hand side of the equations (2.3) for 
a general (/j-dependent metric, gab = 9abit, p, (f), and solve the equations in terms which are 
~ 0{e'^). However, we are primarily interested in the rotation of inertial frames at the axis 
where uj is not (/^-dependent. Therefore we calculate the average of uj over if at each radius. 
To find this we consider the average over (p of the right-hand side of (2.3) and obtain an 
effective matter source that is axially symmetric of the form 

f-2n 

<Sab =< dai^dbip >= / dai^dbixiip. (5.3) 

Jo 

Since the source becomes then axisymmetric, we consider the metric gab^ resp. ^ab axisym- 

metric as well. However, the Killing vector T} is not in general hypersurface orthogonal, and 

we cannot use the vacuum equations to simplify the choice of coordinates. A sufficiently 

general axisymmetric 3-metric has the form 

da"^ = e^^ [dt^ - dp^) - W'^ {d^ - ivdtf , (5.4) 

where 'j,W and to are functions of t and p. Notice that this is precisely the 4-metric (4.7) 
reduced to the 3 dimensions and rescaled by the norm of the Killing vector ^. The metric 
(5.4) satisfies the relations 500 = "511) 501 = 512 = 0, which can be achieved by a suitable 
choice of coordinates. The Ricci tensor components for the exact metric (5.4) are listed in 
Appendix A. Within our approximation scheme it is sufficient to consider only terms ~ O^e^). 
Since for e = the metric (5.4) becomes flat, we can write 

W{t,p)=p + ^w{t,p), (5.5) 
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where w ~ 0{p^) at p ^ 0, but it will be convenient to work with W . Henceforth, e^ is 
absorbed in the corresponding symbols, so it is assumed that -;/; is of the order ©(e), and 
7, w and oj and their derivatives are all of C'(e^); their products can be neglected. From the 
expressions for TZab in Appendix we get, to order ©(e^), the following set of equations: 

7^oo = -i + i' + -i --W , nm = -i--w' , 7^o2 = ;r^(pV)', (5.6) 

P P P P 2p 

7^ll = J - -f" + -j' - -w" , 7^l2 = ^pV , TZ22 = p{w-w"). (5.7) 
p p 

The wave equation (2.4) for the field tp(t, p, cp) becomes, in our approximation, the flat space 
wave equation in cylindrical coordinates: 

i^-r- -^' - \dl^ = 0. (5.8) 

P P^ ^ 

For W = p, cj = we recover the Einstein-Rosen waves. In their case one can determine 
7(t, p) by a quadrature. Now the equations 

TZab = 2 < dai'dhtP >= 2Sab, (5.9) 

with TZab given by (5.6) and (5.7) should determine the functions 7, W and to. 

It is well-known that the gravitational field in the 3-dimensional spacetime has no dynam- 
ical degrees of freedom; these are all contained in the matter fields, in our case in the field 
V'. The entire curvature (5.1) is determined by the local distribution of tp. Outside matter 
the spacetime is flat. In the canonical gravity language, the 2-dimensional metric and the 
conjugate momenta have six independent components, three of which can be annulled by the 
choice of a spacelike hypersurface (time) and two suitable coordinates on the hypersurface. In 
addition, there are three constrained field equations (see e.g. [39] for an interesting exposition 
of Einstein's theory in a 3-dimensional spacetime). 

We have chosen suitable coordinates in the metric (5.4) already. Hence, we have to look 
at the constraints. One is evident. The equation 

^02 = Y (/^^^')' = 2 < i^d^i^ >, (5.10) 

following from (5.6), gives with suitable boundary conditions the 'angular velocity' uj{t,p) in 
terms of the averaged 'angular momentum' component S02 of rotating ^- waves. This equation 
will be analyzed in depth in the accompanying paper [14]. In the explicit construction we 
shall demonstrate how rotating gravitational waves rotate local inertial frames in a 'global 
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sense', like rotating 'ordinary' matter does. Alternatively, the integration function of time is 
fixed by the equation TZi2 = 25i2, see (5.7). The equation for uj still leaves to undetermined up 
to an arbitrary function of time uJo{t). This arbitrariness, however, corresponds just to going 
over to rotating axes (cf. the end of Section 4). Another constraint equation, TZqi = 25oi, 
see (5.6), can be written as 

-J^ = 25oi where J=- = j-W'. (5.11) 

P 

This becomes (2.8) for Einstein- Rosen waves (W = p). 

There is still a constraint equation, the one corresponding to the Hamiltonian constraint 
in canonical gravity. Its left-hand side is given by T^oo ~ 29^^'^ which, in our approximation, 
is equal to 2(^00 + ^11 + P~'^'^22)- Substituting for the Ricci tensor from (5.6) and (5.7), 
we obtain the constraint in the form 

-^' = -(7 - Wy = Soo + Su + ^522. (5.12) 

P P P^ 

Before we indicate how to solve the field equations, we investigate the consistency of 
(5.10) and (5.11). This should generalize the condition (7)' = (7')' for the Einstein-Rosen 
waves when 7 is expressed from (2.8), whereas 7' from the sum of (2.6) and (2.7). There the 
consistency is guaranteed by the wave equation (2.9) for the source V- Now ip satisfies the 
flat-space wave equation (5.8) in cylindrical coordinates with ^^p^ 7^ 0. As a consequence of 
this equation the effective energy-momentum tensor 

Tah = da'4)db%lj-\gagj'"^dcipdd'4) , 5^5 = (1, -1, -p^), (5.13) 

satisfies the conservation law 

VbT^ = -^db (y^T^) - iT-'^dbg.d = 0. (5.14) 

Denoting by Sab = datpdbiJ without the averaging so that < Sat >= Sab, we easily get the 
relations 

Tab = Sab - Isgab => Sab = Tab-Tgab where s = -2T = V'^-^'^- ^(a<pV)^- (5.15) 
The conservation law (5.14) written in terms of Sab reads explicitly as follows 

Soo soi - s'qi ^9^so2 - U = 0, (5.16) 

p p^ ^ 

soi sii - 4i 7^d^si2 + -^S22 - hs' = 0, (5-17) 

p p^ p-^ '^ 

S02 -512 - S12 2^ifS22 - \difS = 0. (5.18) 
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The averaging over (/? makes the dipSab-teTins to vanish so that Sab satisfies the equations 

^00 



-Soi — Sqi - 


-',^ - 


= , Sqi Su - 

P 


- 5n + ^522 - 


-¥' 


= 0, 


(5.19) 


Sqi 5i2 - 

P 


" '-'12 - 


= with S = Soo - 


- Su 5-522- 

p^ 






(5.20) 



The compatibihty of (5.11) and (5.12) requires 

[T)' = 2{pSoi)' = [T'y = p{Soo + Su + ^522). (5.21) 

This is indeed satisfied as a consequence of (5.19). 

What boundary conditions do we impose? As emphasized above, all the dynamics is gov- 
erned by the field ip^t, p, (p) satisfying the -wave equation (5.8). We consider decaying solutions 
at spatial infinity -which guarantee that the 3-dimensional spacetime is asymptotically flat 
although the metric at spatial infinity has a "conical singularity" "which is the measure of the 
total energy of the scalar field ijj computed using the Minko-wski metric. This -will generalize 
naturally the case of the Einstein- Rosen -waves (see [15] for the details). Hence -we assume 
the solutions of (5.8) to have the expansion on each Cauchy surface t=const of the form 



V'(t,p,v?) 



p 



1/2 



/o(t,(p)+X; 



fk{t,<f) 



pk 
k=l ^ 



(5.22) 



The solutions of the -wave equation in 2-|-l dimensions -well-behaved at infinity involve the 

Bessel functions -which decay as p^^'^. Our explicit solutions constructed in the accompanying 

paper [14] -will all admit the expansion of the form (5.22). 

Since the spacetime should be asymptotically flat at infinity the function W must behave 

as follo-ws: 

Wit,p) = p + woit) + f2'^. (5.23) 

fc=i ^ 

However, regarding the field equation (5.7), 

^22 = p(W - W") = 2522, (5.24) 

we see that wo{t) = wi{t) = since S22/P ~ p"^ as a consequence of (5.22). Hence, we may 
write 

Wit,p) = p + f2'^. (5.25) 

fc=2 ^ 
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A similar argument based on the field equation 7^o2 = 25o2 implies the following behavior of 
the 'dragging factor' uj: 



w(t,p) = — ^ + 2^-fe^- (5.26) 

^ k=l ' 

Since S12 ~ p^^, the equation TZi2 = 25i2 with TZi2 given in (5.7) implies 

Woo(i) = const. (5-27) 

This constant is directly related to the total angular momentum in the waves [see equa- 
tion (4.9)]. 

Now in order to guarantee a regular axis, without a conical singularity present, we require 
the following conditions at p = 0: 

W(t,0) = p , VF'(t,0) = l , u){t,Q) = ujQ{t) , u;'(t,0) = , 7(t,0) = 0. (5.28) 

The only missing boundary condition is that on 7 at infinity. To find it, we return back to 
the "Hamiltonian" constraint, multiply by p and integrate to obtain 

-i{t,p) -^{t,d) -W'{t,p) + W'{t,d) = 1% ( Sm + 811 + ^822] dp + \{t), (5.29) 

where \{t) is an arbitrary function of time, and 7(t, 0) = 0, VF'(t,0) = 1 by the boundary 
conditions (5.28). Regarding the behavior of W{t,p) in accordance with (5.25) and realizing 
that Sqi =< VV'' >~ P~^) ■we can deduce from the constraint (5.11) that \{t) must be a 
constant. Since, in addition, we have to satisfy the boundary condition 7(t, 0) = 0, this 
constant has to vanish and the integrated constraint equation (5.29) becomes 

7(t, p) - W'{t, p) + l= Tp (soo + 5ii + 4^22 j dp. (5.30) 

At p = 0, W'{t, 0) = 1 and, indeed, -/{t, 0) = 0. Now for p ^ 00, W'{t, cx)) = 1 + 0{p~^), [cf. 
(5.25)], so putting 7(t, 00) = 700, we get 

7oo = ^" (500 + 5n + ^522) pdp = j^ l^<^^> + < ^'^ > + <(^^V')^> ^ ^^^^ 

(5.31) 
Therefore the function 7 at spatial infinity approaches a constant which is equal to the total 
energy of the scalar field tp computed using the Minkowski metric. In the case of the Einstein- 
Rosen waves, the term p^^ < (d^ip)'^ > vanishes, and we recover the result (2.35) in [15]. 
Although 7oo is energy for weak fields, the physical Hamiltonian turns out to be a non- 
polynomial function of 700 for in general strong fields [40], [15]. 
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In order to find the complete solution of the averaged Einstein equations up to order 

0{e'^) we have to solve the equation (5.24) with (5.28), and with the right-hand side given in 

terms of ip. The function 7 can then be determined from (5.30). Alternatively, we can take 

the combination ^(— T^oo + "^11 — P^^T^22) and obtain an equation for 7, similar to (5.24) for 

W, in the form 

7 - 7" = cSii- 500 + 4^22. (5.32) 

/?^ 

Both equations (5.24) and (5.32) are one dimensional wave equations with given right-hand 
side and boundary conditions 7 = 700, W = 1 at p ^ 00 and 7 = 0, VF = 1 at p = 0. Both 
equations have thus the same form as the equation for the motion of an elastic string with 
fixed end points under the influence of an external force. There are standard techniques for 
solving such equations. For example, one can find the desired result by expanding the solution 
in time t in terms of eigenfunctions sin np and expanding the "external force" similarly (see 

[41]). 

To determine the "dragging potential" w(t, p) we employ simpler equations 7^o2 = 25o2 
and 7^12 = 25i2, with TZq2 and 7^12 given by (5.6) and (5.7). The explicit solution of 
these equations with the source field ip representing specific ingoing and outgoing pulses of 
the rotating cylindrical gravitational waves demonstrates nicely the dragging effects of such 
waves. We treat this problem in-depth in paper [14]. 

6 Dragging of inertial frames in axisymmetric spacetimes with 
cylindrical waves and spinning string 

The vacuum metrics (4.7) considered by Mashhoon et al [31] as "rotating gravitational waves" 
cannot have a regular axis since this requires the Killing orbits to admit orthogonal surfaces. 
However, the axis can represent a rotating cosmic string. Let us assume the 'rotation pa- 
rameter' uj in (4.7) is small so that the terms in 0{uP') can be neglected. The inspection 
of vacuum field equations following from the Anzatz (4.7) then reveals that one can choose 
W = p and equations for Tp and 7 are the same as (2.9) and (2.10), see Appendix A. The 
rotational perturbation uj is determined by the equation 

i?12 = -^(A-2V)-=0, (6.1) 



and the constraint equation 



i?02 = -^(A-2V)' = 0. (6.2) 
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These two equations are directly related to the properties of the angular momentum of 
cylindrical systems derived from the Komar expressions (4.1) or (4.3) and (4.6). In (4.6) 
we indicated that one can evaluate the integral at p ^ oo since its value in vacuum does not 
depend on the location of the cylinder C. For the metrics of the form (4.7) with small uj we 
obtained the result (4.9) which again one can evaluate at p ^ oo. However, it is evident from 
the constraint equation (6.2) that the same value for J results at any p ^ since p'^e^'^'^uj' is 
independent of p and by (6.1) independent of time. 
Integration of (6.1) and (6.2) yields 

g27(t,p) ^ 

to' = 5 J, (6.3) 

p-^ IT 

where the integration constant was chosen in accordance with the expression (4.9) for the 
angular momentum. The angular momentum is a small constant parameter with — having 
the dimensions of length. Hence, 

uj{t,p) = J ^dp + iVo{t); (6.4) 

we put ujQ = because it can be transformed away by going to rotating axes, without a 
change of J (cf. Section 4). 

When no wave is present, (6.4) gives oj = los = ~^^j the metric (4.7) with W = p 
becomes precisely the metric (4.11) of the massless spinning string with angular momentum J. 
The "string perturbation" ujs causes rotation of local inertial frames at some finite p as 
compared with frames at infinity. However, this relative rotation is time-independent as, 
e.g., in the Kerr metric. 

Now consider -0 and 7 to represent an exact cylindrical gravitational wave, e.g., the 
Bonnor/Weber-Wheeler wave given in (2.12) and (2.13). Since the rotational effects are 
given by a linear perturbation of the metric we can write the total perturbation (6.3) as 
uj' = Lu'g + uj[y, the time-independent uj'g is attributed to the spinning string, whereas the 
time-dependent part 

-w = ^(e^^(*-) - 1) (6.5) 

is associated with the presence of a cylindrical gravitational wave interacting with the string. 

Now at t"^ ^ a^ + p^, (2.13) implies 7 ~ 0, so when the pulse is far away from the axis, 

the dragging near the axis is dominated by the string. As the wave proceeds towards the 
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axis, the dragging due to its presence increases. The expansion of 7 at small p [see (2.11)] 
reads 



J2 



t 



7 = 86V^^^^ + 0(p^). (6.6) 



^'w-^-^^^TTT^^A^ (6-7) 



^w - _„^ In p-rrr^^- (6-8) 



1 

7 



1 - ^ + 0(p 



^4\ 



(6.9) 



Hence, at small p equation (6.5) gives 

1662kJ 1 P 
7ra3 ^(1 + ^2)4^ 

so that 

1662kJ, ^ P 

Since we assumed uj - or, correspondingly, J - small, our approximation will break down^ at 
p = 0, but it should be reliable at all p's for which Kj/irp ^ 1; this can be assured at any 
p 7^ by choosing J sufficiently small. 

At spatial infinity, t =const., p -^ 00, the function 7 approaches a constant 

2 
h^ 1 - 
2 

the time-dependent terms appear only in 0{p~'^). Substituting this expansion into (6.5) and 
integrating, we find the angular velocity of the rotation of inertial frames at large p: 

^w = ^^(e^'-l)+0{p-'). (6.10) 

The perturbation los corresponding to the string also decays at infinity, so local inertial 
frames do not rotate with respect to the 'fixed stars' at infinity, i.e., with respect to the 
'lines' ip =const. However, they do rotate with respect to these lines close to the axis, at 
small /?, with the time-dependent angular velocity oow given by (6.8). This vanishes at times 
|t| S> 1 when the incoming or outgoing pulse is far away from the axis. It becomes maximal 
{uj\\/ = 0) at |t| = l/\/3. Then the pulse is close to the axis. 
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*A somewhat analogous situation arises in the problem of slowly rotating collapsing spherical shells when 
the shells approach the horizon [8] . 
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A The Ricci tensor for the averaged rotating gravitational waves 
in 2+1 dimensional spacetime 

Assume the metric is of the form (5.4). The Ricci tensor components read as follows: 

„ W W , W 

^00 = -7 + 7 -^ + ^1+^1 

+ e-^^W h-f'WiOLo' + Wlo^ - W"io'^ - SW'ujuj' - \WJ^ - WujJ' - \W^uj'^J^e~'^A , 

W W' W , n / 1 o • \ 

7^0l = —'r-^ + i^7+e-^^W[W^iOLo'-lWLOi:o'-^WLOLo'j, 

7^o2 = e-^^/wdWio' + lWio" + W"io-W-f'uj' -Wio + lW^iOio'^e-^A , 



nn = 7-7" + ^7'-^ + ^7 + iW^VV^ 



7^l2 = y-^^'W (Wu' + 3Wuj' - 2Wjuj'] , 
7^22 = e-^^w(w-W" -IW^co'^e-^A. 



2 ' 
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